Level density of a Fermi gas and integer partitions: a Gumbel-like finite-size correction 
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We investigate the many-body level density of gas of non-interacting fermions. We determine its 
behavior as a function of the temperature and the number of particles. As the temperature increases, 
and beyond the usual Sommerfeld expansion that describes the degenerate gas behavior, corrections 
due to a finite number of particles lead to Gumbel-like contributions. We discuss connections with 
the partition problem in number theory, extreme value statistics as well as differences with respect 
to the Bose gas. 
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I. INTRODUCTION 

The many-body (MB) level density is a fundamen- 
tal quantity which describes basic properties of thermo- 
dynamical systems. Its behavior is of principal interest 
in many different fields such as nuclear decay rate [1, 2], 
thermonuclear reactions in stellar processes or - more re- 
cently - black hole entropy [3-5] . Unfortunately ab initio 
computations are out of reach for interacting systems and 
for large particle numbers. However a statistical treat- 
ment is still possible when interactions are taken into 
account in a mean field approximation. For fermions, 
two extreme regimes can be identified. At temperatures 
which are low compared to the Fermi energy /ip but high 
with respect to the mean level spacing at no (degener- 
ate gas approximation), has a stretched exponential 
growth [6] 
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where Q is the excitation energy and p is the smooth 
single-particle (SP) level density. Nowadays Bethe's the- 
ory still remains a basic ingredient for advanced MB 
theories like back-shifted Fermi gas models [7] and the 
Gilbert-Cameron model [8]. From Eq.(l), we see 
is essentially sensitive to the SP level density around 
the Fermi energy. The second regime corresponds to 
high temperatures Q S> /ioj where the classical limit 
is reached and where has a power-law behavior 

(Maxwell-Boltzmann regime). The aim of the article is 
to analyze, for fermions, the transition between the two 
regimes (see Fig.l). To tackle this, we compute - in the 
degenerate gas limit - the effect of a finite number of par- 
ticles when the temperature increases. For a power-law 
SP level density p(/io), we find two contributions to p^^g: 
the first one is given by the usual Sommerfeld expansion 
while the second one is related to low-lying level exci- 
tations. This work is in connection with similar results 
obtained by Comtet et al for bosons [9] . In a mean-field 
approximation, p^^ is also directly related to a combi- 
natorial problem: to find the number of configurations 
of the SP occupied levels for a given excitation energy 



[10, 11]. 

The paper is organized as follows. In Sec. II we introduce 
basic definitions of the MB level density. We calculate the 
expression of the smooth part of the MB level density in 
the low- and high-temperature regimes. We establish a 
relation between the partition problem in number theory 
and the MB level density by using an equidistant spec- 
trum in Sec. III. We compute corrections due to a finite 
number of particles from an exact asymptotic formula. 
Finally we generalize in Sec. IV the result to a power-law 
spectrum using a non-local Sommerfeld expansion pro- 
vided by Garoni et al. [12]. 



II. GENERAL FRAMEWORK 

The MB level density is a discrete quantity which, for 
a gas in a SP potential, is given by: 
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It is a microcanonical quantity where the particle number 
N and the energy of the system E are fixed (all quanti- 
ties are dimensionless). Cj corresponds to the energy of 
the jth level of the SP potential for one particle and mj 
is the occupation number of the jth SP level for a given 
configuration {rrij}. In the case of fermions without spin 
degeneracy, the Pauli principle leads to nij being either 
or 1 . We define Q as the difference E — Eq, where Eq 
is the ground state energy. For each configuration {mj}, 
the conservation of particle number and energy yields 
J^JLo ™j ~ ^ ^i^d J^JLo ^j'^j ~ ^1 respectively. 
Now we write Eq.(2) in term of an inverse Laplace trans- 
form: 
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where 



5(/3,p) = /3 
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is the entropy. The constant b (resp. c) is chosen such 
that it is larger (resp. smaher) than the sum index 
of exp[— P^l{f3 , ^)] with respect to the variable f3 (resp. 
/i). Both parameters (3 and /i correspond respectively to 
the inverse temperature and to the chemical potential. 
il{/3,IJ.) is the Grand potential defined as 
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where p(e) = X^jlo ^{^ ^ ^j) is the SP level density. The 
saddle point method with respect to and (3 gives the 
following equations: 
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I?(/3, fi) is the determinant of the second derivatives of 
iS(/3, fi) which is not relevant for the following discus- 
sion. The last two equations define ij,{E, N) and /3{E, N). 
From Eq.(7), we see that fi must decrease from positive 
to negative values when the temperature increases. Ig- 
noring the discreteness of the SP level density, we now 
replace p by a smooth function: 



p{e) « p{e) = ve"' 
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For instance, the case v — D \s the leading order of the 
i^-dimensional HO. The case of billiards is given by = 
D/2, where D is the dimension of space. When p is given 
by Eq.(9), equations Eqs.(7) and (8) become: 
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where Li,y(x) = 'Yl'k^i^^ /^^ ^he polylogarithm func- 
tion. In this case the total energy and the Grand poten- 
tial are connected by i? = —M{l3,p). From Eq.(4), the 
entropy simplifies into 



S = /3 



{l + l/iy)E - pN 



(12) 



For the low temperature regime, we expand Eqs.(7) and 
(8) with l/{Pp) as a small parameter. This so-called 
Sommerfeld expansion [13, 14] is valid for any regular 
fimction /: 



(13) 



where Q = T]{2j)d^^^~'^^ /(e)/de(2i-i) |^^^ and 7](x) is 
the Dirichlct eta function. Using Eqs.(7), (8) and (13), 



we define the chemical potential and inverse temperature 
to leading order - pq and /3q respectively - as: 
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: E-Eo = 27K2)p(mo)/3o" , 

where Eq = ep{e)de. We have neglected terms in the 
derivatives of p in Eqs.(14) and (15). When p satisfies 
Eq.(9), we can easily invert (14) and (15) to obtain: 
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and 



/3o = z^V2^(2)/QiV^""'^/'^ 

(16) 

Integrating the Grand potential Eq.(5) by parts, we get: 
nipo, /3o) = ^^(0, Mo) - 2f]{2)p{po)l3o^. (17) 

Using the expansions of Eqs.(14), (15) and (17) in the 
entropy Eq.(4), the result is 



5(/3o,/^o) = /So 



'n{Q,po) + 2r^{2)p{po)l3-^ 



At zero temperature it leads to: 

-n{0,po)+Eo- poN 
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-E-poN . 
(18) 
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We have used the latter equation in Eq.(18). Finally we 
obtain the leading order of the entropy for a Fermi gas 
at low temperature: 



S{l3o,po) = 2\j '—p{po)Q , 



(20) 



where we have used ri{2) = 7r^/12. This result holds 
for any smooth SP level density; in the case Eq.(9) 
the entropy is written S{Q,N) = 2^iyTT'^Q/6N'-''-^'^/^'' . 
Higher order terms can be easily computed leading to 
a power series in Q and N. We emphasize that for 
fermions, the dependence of S on the excitation energy 
docs not depend on the system. In contrast, its N de- 
pendence comes through p{pq). This behavior is very 
different from the bosonic case (except for ly — 1, see be- 
low) where the energy variation of the entropy strongly 
depends on the SP level density, i.e. p^g oc exp[(3'^/*^'^"'"^^] 
[9, 10]. 

When the particle number is fixed and when the temper- 
ature is sufficiently large, Eq.(7) leads to a large negative 
chemical potential, e^^ tends to and the Fermi-Dirac 
statistic converges to the Maxwell-Boltzmann distribu- 
tion. In this case Eqs.(lO) and (11) become: 
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We invert Eqs.(21) and (22) to get p and /3 as a func- 
tion of and E. Using Eqs.(6) and (12), it yields the 
following formulas for E' ^ A^: 
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and = liilP" N /V{v+\)\/ P . We also find tlie usual equa- 
tion of state for perfect gases: E = vNT. In contrast to 
the stretched exponential growth at low temperature, we 
notice that the MB level density has a power-law behav- 
ior at high temperatures. 



III. THE ONE DIMENSIONAL HARMONIC 
OSCILLATOR 

In this section we compute exactly the MB level den- 
sity for a system of non-interacting particles confined by 
a one-dimensional harmonic potential (IDHO) with unit 
frequency and make connections with number theory. We 
choose £0 = such that the spectrum is given by integers. 
We must compute the MB level density for a given inte- 
ger excitation energy n and for a given particle number 
N , by counting all possible configurations of particles. 
A schematic way to depict this enumeration consists of 
using Young diagrams. They form n boxes with j non- 
increasing lines such as: ci > C2 > ... > Cj , X]i=i ~ n 
where a is the column number for the line i with j < N . 
Consequently, for an integer n the number of Young dia- 
grams is the number of ways of counting n as non-zero in- 
tegers when the line number is less than N. This number 

(n) is the so-called partition function of n restricted to 
at most N terms. 

For each Young diagram we associate a unique config- 
uration of particles in such a way that the horizontal 
axis represents the excitation energy of the particle and 
the vertical axis represents the index of the occupied 
level. Hence we have a one-to-one correspondence be- 
tween each particle configuration and Young diagrams; 
Pfi{n) is p^g{n,N) [10, 11]. In the same way one can 
easily shows that the result holds for particles with Bose- 
Einstein statistics as well. 

The theory of partitions is studied in profusion by math- 
ematicians (see [15] for an introduction). One of the 
most singular results comes from the works of Hardy, Ra- 
manujan and Rademacher who computed an analytical 
expression when the number of terms into p„ (n) is not 
restricted anymore [16, 17]. In this ca.se p^^^[n) = p{n) 
is called the partition number of an integer n. The lead- 
ing asymptotic expression of p{n) for large n is given by 
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To leading order the exponential term in Eq.(24) is equiv- 
alent to Eq.(6) with Eq.(20) and = 1 as expected. The 
correction term for a finite number of terms to p„ (n) was 
computed by Erdos and Lehner [18]: 
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with.g = ^7r2 /{6n)N + l /2 ln[7rV(6ri)] and iV = o{n^/^). 
This distribution is the so-called Gumbcl-likc, and is one 
of the universal distributions in the domain of extreme 



values statistics [19]. It shows how the finite particle 
number contributes to the decrease of the number of con- 
figurations when heating the system (increases n) and 
determines the transition from low- to high-temperature 
regime. The main target of the article is to evaluate this 
correction when the spectrum is no longer equidistant 
but fulfills the more general law Eq.(9). This will be 
computed in the following section. 

Note that in the case of fermions, the Fermi energy 
depends on the particle number, while this is not the 
case for bosons. This fundamental difference implies that 
the equivalence between the two statistics is only valid 
in the case of the IDHO. In the context of the partition 
theory the generalization to a SP level density of the form 
Eq.(9) leads to map p^g{E, N) to the partition number 
of an integer n = E into a sum of N distinct (due to 
Pauli principle) 1/v-th, powers of integers. For bosons, a 
similar relation exits but without any distinction between 
summands [9, 10]. 



IV. GENERAL CASE 

In the degenerate gas limit the chemical potential is 
large with respect to the excitation energy. At low tem- 
perature, the fugacity z = e^^^ tends to infinity. To solve 
the equations Eqs.(lG) and (11), we must find the asymp- 
totic expansion of the polylogarithm function for large 
negative argument. Garoni et al. give the inversion ex- 
pression of Li^(—z) [12]: 

Li.(-z) = -2 "^^^"^/"T'"!!"'^ -cos(7rz.)Li.(-z-i) 



ri=0 



r(i^-|- 1 - 2n) 



(26) 

Note that for v half-integer, the term in Li^{~z ^) van- 
ishes. Furthermore when v is an integer, the sum in 
Eq.(26) tends to the integer part of i'/2. At low temper- 
ature tends to zero and limLi,y(.T) = x. Eqs.(lO) 

and (11) with Eq.(26) become: 



^ = 2r(i. + i)^ 



— 271 ,,iy — 2n 



r;(2n)/3-^"^ 
r(i^ + 1 - 2n) 



-cos(7r.) ^^"+j,^ (27) 
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and 



E = 2 iyT{iy+l)Y^ 



77(2n)^-^"// 



2n,,i/+l-2n 



n=0 



r{iy + 2- 2n) 



+ cos(7rt.) ^ . (28) 



Terms in the sum on the r.h.s. of Eqs.(27) and (28) corre- 
spond to the Sommerfeld expansion terms. The meaning 
of the exponential term in Eqs.(27) and (28) will emerge 
from the following example. With = 1, Eqs.(27) and 
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FIG. 1: Upper panel: The MB level density for i/ = 1 and for 5 
fermions as a function of the excitation energy Q. The dotted 
line is the exact computation. The full line corresponds to the 
degenerate or low temperature regime, Eq. (24) (with n = Q). 
The dashed line corresponds to the high temperature regime, 
Eq. (23) (with E ~ Q; we have included the expression of the 
determinant T> = [see Eq.(6) and below for discussion]). 
Lower panel: The function F{Q, N) for the same system. The 
dotted line corresponds to the exact computation and the full 
line to the analytical expression Eq.(34). 



(28) yield: 



N 
E 



/3 ' 
2 /32 



/32 



(29) 
(30) 



The computation of S to the second order is done by 
expanding /i and (5 in Eqs.(29) and (30). We have: 
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-^2r,(2)/QW 



(31) 



(32) 



Using Eqs.(29), (30), (31) and (32) into Eq.(12), we get 

-t 



S = 2v/2r/(2)Q 



(33) 



where t = ^J2'q{2)/QN + h\[yj2T]{2) / Q]. Subtracting 
the leading order entropy Eq.(20) to Eq.(33) and using 
Eqs.(6) and (9), Eq.(25) is found with n^Q. The Gum- 
bel law is thus recovered, as already stated in the pre- 
vious section. The main contributions to are given 
by the Sommerfeld terms. These terms depend only on 
the local properties of the SP level density close to the 
Fermi energy, and ignore the fact that the single parti- 
cle spectrum is bounded from below. The transition to 
the MaxwcU-Boltzmann regime, where all particles con- 
tribute to a single configuration, cannot be described by 



this expansion. The exponentially small term e~^^ in 
Eqs.(27) and (28) takes into account this effect, and this 
is the reason why we call it a non-local correction to the 
MB level density. For a given z^, we compute the non- 
local contribution to p^^ using the previous method, to 
obtain 



F{Q,N) ^ ^ 



exp 



cos(7rz^)r(i' 



■ exp(-/3oMo) 



(34) 



JJ^IQ, N) represents the smooth MB level density com- 
puted by means of the (local) Sommerfeld expansion 
Eq.(13). Equation (34) is the main result of the pa- 
per. The final correction takes the form of a modified 
Gumbel law. This is different with respect to the case 
of bosons, where, depending on the value of j/, three 
different functional dependencies of the correction were 
found [9]. They correspond to the three different dis- 
tributions that appear in the theory of extreme value 
statistics. The difference is related, in particular, to the 
fact that there is no Bose-Einstein condensation for non- 
interacting fermions. 

The qualitative behavior of F{Q,N) has a strong de- 
pendence on the value of ly. For instance, F{Q,N) = 1 
for V half-integers. It can be larger than one for some 
specific values of v. This effect is not surprising, since 
for ly =^ I the Sommerfeld expansion gives the leading 
variation in N. In that case F{Q, N) is not constraint 
to be less than one. The behavior of the MB level den- 
sity as a function of the excitation energy is illustrated 
in Fig. 1 for i/ = 1 and for TV = 5 fermions. This value 
of v describes a one-dimcnsional harmonic oscillator, but 
also for instance a two-dimensional billiard of arbitrary 
shape. A low number of particles has been chosen in or- 
der to enhance the finite N cotrrections. The upper panel 
illustrates the transition from the degenerate gas regime 
to the Maxwell-Boltzmann one. The lower panel shows 
the correction factor F{Q,N). We note that the correc- 
tion term is significant only in a small region around the 
degenerate regime, i.e. for Q ^ 15. At higher excitation 
energies, the Maxwell-Boltzmann regime is reached. 

In conclusion we have demonstrated that in the tran- 
sition from the low to the high temperature regimes the 
MB level density of a Fermi gas is described by a mod- 
ified law of extreme value statistics. This result is valid 
for a large set of systems according to Eq.(9). In the 
context of number theory, our results apply in particular 
to the partition of integers as a sum of powers of distinct 
integers. Our results are of interest in the computation 
of the level density of small nuclei for which the finite 
size correction term (34) is the most important. 
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